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Abstract 

In this article, we construct the diquark-diquark-antiquark type interpolating currents, and 
study the masses and pole residues of the = | and I”*" hidden-charm pentaquark states 
in details with the QCD sum rules by calculating the contributions of the vacuum condensates 
up to dimension-10 in the operator product expansion. In calculations, we use the formula 
fi = ^Mp^ — (2Mc)^ to determine the energy scales of the QCD spectral densities. The present 

predictions favor assigning the Pc (4380) and Pc (4450) to be the | and pentaquark states, 
respectively. 
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1 Introduction 

In 1964, Gell-Mann suggested that multiquark states beyond the minimal quark contents qq and 
qqq maybe exist [1] , a quantitative model for the tetraquark states with the quark contents qqqq was 
developed by Jaffe using the MIT bag model in 1976 [5]. Latter, the five-quark baryons with the 
quark contents qqqqq were developed [5] , while the name pentaquark was introduced by Lipkin [3] . 
The QCD allows the existence of multiquark states and hybrid states which contain not only quarks 
but also gluonic degrees of freedom. We can construct the tetraquark states and pentaquark states 
according to the diquark-antidiquark model and diquark-diquark-antiquark model, respectively 
[5l|6]. In the light quark sector, the nature of the scalar mesons below I GeV is under controversy 
[7], although those light tetraquark states are not ruled out in the limit [5]. In the heavy 
quark sector, several X, Y and Z mesons are observed, such as the ^^(3900)^, ^^(4020/4025)^, 
Z(4430)^, the net charge indicates that their constituents are ccud or ccdu, for recent review on 
both the experimental and theoretical aspects, one can consult Ref. [5]. Some X, Y and Z mesons 
are assigned tentatively to be tetraquark states, irrespective of the diquark-antidiquark type or 
the meson-meson type. The two heavy quarks play an important role in stabilizing the multiquark 
systems, just as in the case of the (^“e+)(/r+e“) molecule in QED [TU]. The spacial separation 
between the diquark and antidiquark in the tetraquark states [lOl HI) (or meson and meson in 
the molecular states [HIS]) may lead to small decay widths, we can study the decay patterns 
by performing the Fierz rearrangements non-relativistically in the Pauli-spinor pace [n Hairs] or 
relativistically in the Dirac-spinor space [T4] . 

Recently, the LHCb collaboration observed two exotic structures (A: (4380) and Pc(4450)) in 
the J/ipp mass spectrum in the A® —^ J/ipK~p decays, which are referred to be charmonium- 
pentaquark states now [TS]. The Pc(4380) has a mass of 4380 ± 8 ± 29MeV and a width of 
205 ± 18 ± 86MeV, while the Pc(4450) has a mass of 4449.8 ± 1.7 ± 2.5 MeV and a width of 
39 ± 5 ± 19 MeV. The preferred spin-parity assignments of the Pc(4380) and Pc(4450) are = 
I and 1^, respectively. The significance of each of the two resonances is more than 9 ct [T5] . 
The Pc(4380) and Pc(4450) have attracted much attentions of the theoretical physicists, several 
attempted assignments are suggested, such as the YcD*, E*D*, xaP molecular pentaquark states 
m (or not the molecular pentaquark states m), the diquark-diquark-antiquark type pentaquark 
states [18] , the diquark-triquark type pentaquark states m, re-scattering effects [20], etc. We can 
test their resonant nature by using photoproduction off a proton target [21j . 
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The quarks have color SU (3) symmetry, we can construct the pentaquark states according to 
the routine quark —>■ diquark —> pentaquark, 

(3 (g) 3) 0 (3 0 3) 0 3 = (3©6) 0 (3© 6) 0 3 = 30 303© •• • = 1 © • •• , (1) 

or construct the molecular pentaquark states according to the routine quark —>■ meson and baryon 
—>■ molecular pentaquark state, 

(303)0(30303) = (1© 8) 0 (1 ©•••) = (1 ® 1) ©••• = 1 ©••• , (2) 

where the 1, 3 (3), 6 and 8 denote the color singlet, triplet (antitriplet), sextet and octet, respec¬ 
tively. In the diquark model, the pentaquark states consist of two diquarks and an antiquark, 
which are colored constituents, it is easy to form compact bound states due to the strong attrac¬ 
tions at long distance. In the meson-baryon model, the molecular pentaquark states consist of 
a colorless meson and a colorless baryon, attractions induced by exchanges of the intermediate 
mesons (Yukawa-like potentials) are needed to form loose bound states. In this article, we take 
the Pc(4380) and Pc(4450) as the diquark-diquark-antiquark type pentaquark states, construct the 
interpolating currents consist of five quarks according to Eq.(l), and study their masses and pole 
residues with the QCD sum rules. 

In previous works, we described the hidden charm (or bottom) four-quark systems qq'QQ by a 
double-well potential [Hums]- In the four-quark system qq'QQ, the Q-quark serves as a static 
well potential and combines with the light quark q to form a heavy diquark 2?® g in color antitriplet 


q + Q , (3) 

or combines with the light antiquark q' to form a heavy meson in color singlet (meson-like state in 
color octet) [HJ [53] 


q' + Q ^ g-'Q(g'A“Q), (4) 

the Q"C[uark serves as another static well potential and combines with the light antiquark q' to 
form a heavy antidiquark ©’g^g in color triplet |14] . 

q' + Q T’f,g, (5) 

or combines with the light quark q to form a heavy meson in color singlet (meson-like state in color 
octet) p2l[23] 


q + Q — t Qq (QX^q) , (6) 

where the i is color index, the A“ is Gell-Mann matrix. Then 

^Q+^'Q —>■ compact tetraquark states , 

q'Q + Qq —loose molecular states , 
q'X°‘Q + QX^q —>■ molecule — like states, (7) 

the two heavy quarks Q and Q stabilize the four-quark systems qq'QQ, just as in the case of the 
(/r“e+)(/r+e“) molecule in QED [TO] . 

The hidden charm (or bottom) five-quark systems qqiq 2 QQ can also be described by a double¬ 
well potential by using the replacement, 

<?l+<?2+0 ^ 
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just like the four-quark systems qq'QQ [Him, where the denotes the heavy triquark in 

color triplet, the q' in the bracket denotes that the 'D^q^q^ is in color antitriplet, just like the q'^. In 
the heavy quark limit, the Q-quark (Q-quark) can be taken as a static well potential, the diquark 
quark q lie in the two wells, respectively. 

The QCD sum rules have been applied extensively to study the hidden-charm (bottom) tetraquark 
states [24], however, the energy scale dependence of the QCD spectral densities is not studied. In 
previous works, we studied the acceptable energy scales of the QCD spectral densities for the hid¬ 
den charm (bottom) tetraquark states and molecular (and molecule-like) states in the QCD sum 
rules in details for the first time [HEHEsjESIEe], and suggested a formula 

- (2Mq)2 , (9) 


to determine the energy scales based on the analysis in Eqs.(3-7), where the X, Y, Z denote the 
four-quark systems, and the Mg denotes the effective heavy quark masses [H[l2l[23]. The energy 
scale formula works well for all the tetraquark states, molecular states and molecule-like states. 

In the non-relativistic quark model, the heavy quarks have finite masses, which quantitatively 
affect the spin-spin interactions between the quarks within one diquark or in two different di¬ 
quarks m- In the QCD sum rules, the net effects of the different dynamics are embodied in 
the effect masses Me and M^, respectively, for example, the Zc(3900) and Z;,(10610) can be 
tentatively assigned to be the = 1“' tetraquark states with the symbolic quark struc¬ 


tures 


[cii]s=o[cc!]s=l-[cu]s=l[cd]s=o 

x/2 


and 


[fcM]s = o[fe^]s = l —[fe«]s = l[Ms=Q 

C2 


respectively, where the subscript 


S denotes the spin, the optimal energy scales of their QCD spectral densities are quite different 
MZc( 3900 ) = 1-5 GeV and fJ.Zt{W6io) = 2.7 GeV [HISS], although they are cousins. While in the 
heavy quark limit mg —>■ oo, we naively expect that the two energy scales fJ-Zc{ 3900 ) ^-^d ^J-Zb{l 06 l 0 ) 
coincide. In this work, we extend the energy scale formula to study the diquark-diquark-antiquark 
type pentaquark states, and try to assign the Pc(4380) and Pc(4450) to be the | and I"*" pen- 
taquark states, respectively. 

The article is arranged as follows: we derive the QCD sum rules for the masses and pole residues 
of the Pc(4380) and Pc(4450) in Sect.2; in Sect.3, we present the numerical results and discussions; 
and Sect.4 is reserved for our conclusions. 


2 QCD sum rules for the Pc(4380) and Pc(4450) 

In the following, we write down the two-point correlation functions n^i^(p) and in tbe 

QCD sum rules, 

n^Ap) = ijAxe^P-A0\T{JAx)M0)}\0), ( 10 ) 

= i J Axe^P-^OlT {J^Ax)Jo.A0)} \0), ( 11 ) 

where 

JAx) = P“e*'^''e'"*"Mj(a:)C75dfe(a;)u^(a:)C'7^c„(x)Cc4(a;), (12) 

JfbAx) = -^e^^‘'e^^^A”^^uJ{x)C-/5dk{x)[uA{x)Cj^Cn{x)-f^Ccl{x)+uA{x)Cj^Cn{x)-ff,Cd^{x)] , 

( 13 ) 

the i, j, k, ■ ■ ■ are color indices, the C is the charge conjugation matrix. The diquarks qJCTq'j, 
have five structures in Dirac spinor space, where CT = Cys, C, Cy^ and for the 

scalar, pseudoscalar, vector, axialvector and tensor diquarks, respectively. The structures Cy^ and 
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are symmetric, while the structures C 75 , C and are antisymmetric. The scattering 

amplitude for one-gluon exchange is proportional to 

= —o {5jk5ii — SikSji) + - {SjkSii + SikSji) , (14) 

where the z, j and k, I are the color indexes of the two quarks in the incoming and outgoing channels 
respectively. The negative sign in front of the antisymmetric antitriplet indicates the interaction 
is attractive while the positive sign in front of the symmetric sextet indicates the interaction is 
repulsive. The attractive interactions of one-gluon exchange favor formation of the diquarks in 
color antitriplet 3c, flavor antitriplet 3/ and spin singlet Is [27] . while the favored configurations 
are the scalar (Cys) and axialvector (Cy^^) diquark states [28ll29]. The calculations based on the 
QCD sum rules indicate that the heavy-light scalar and axialvector diquark states have almost 
degenerate masses [55], while the masses of the light axialvector diquark states lie (150 — 200) MeV 
above that of the light scalar diquark states [29], if they have the same quark constituents. In this 
article, we choose the light scalar diquark and heavy axialvector diquark as basic constituents, and 
construct the scalar-diquark-axialvector-diquark-antiquark type currents (x) and with the 
spin-parity | and respectively to interpolate the pentaquark states Pc (4380) and Pc (4450), 
respectively, see Eq.(3) and Eq.( 8 ). 

In fact, we can also construct the axialvector-diquark-scalar-diquark-antiquark type current 
ri)j.{x) and axialvector-diquark-axialvector-diquark-antiquark type current 







^ila ^ijk ^Imn 

7! 

^ila ^ij k ^Imn 

Te 

+ {fl V) , 


[uJ(x)C'y^Ufc(x)d^(x)C'y 5 C„(x) -I-2uJ(x)C'y^dfc(x)M^(x)C'y5C„(x)] C'cf(x), 

[uj{x)C-ffj,Ukix)d^{x)C-f^Cnix) + 2 uJ{x)Cjf^dkix)u^{x)Cj^c„{x)] ysCcf (x) 

(15) 


to study the spin-parity | and | pentaquark states, respectively. As the masses of the light 
axialvector diquark states lie (150 — 200) MeV above that of the corresponding light scalar diquark 
states [55]. The currents ^^^(x) and rj^vix) are supposed to couple to the pentaquark states with 
larger masses compared to the currents J^{x) and J^v(x), respectively. 

The A° can be well interpolated by the current J(x) = {x)Cj 5 dj{x)bk{x) [30], the u and 

d quarks in the A° form a scalar diquark [ud] in color antitriplet, the decays A° —>■ J/iljpK~ take 
place through the mechanism, 

A°([ud] 6 ) —>■ [ud]ccs ^ [ud]ccuus ^ P^{[ud][uc]c)K~{us) ^ J/ippK~ , (16) 


at the quark level. In the decays P+([u(i][uc]c) —>■ J/ij^p, the scalar diquark [ud] survives in the 

decays, the decays are greatly facilitated. On the other hand, if there exists a light axialvector 

diquark [ud], which has to dissolve to form a scalar diquark [ud], the decays are not facilitated. 

The currents and J^y(O) couple potentially to the ^'*',1 and | , I"*" hidden-charm 

pentaquark states P^, PT and Pt, P 3 , Pt, respectively, 

2 2 2 2 2 


{o\Mo)\Ptip)) 

2 

(0|J^(0)|P3-(p)) 

2 

(0|V(0)|P_+(p)) 

2 

(0|J^.(0)|d^3-(p)) 

2 

(0|d;..(0)|P5+(p)) 

2 


ftp^LU+{p,s), 

2 

s), 

2 ^ 

gtPfiPvU^ip^s) , 

2 

/| [PliPv{P^s)+Pt,U~(p,s)\ , 

A|t/+(p,s), 


(17) 


(18) 
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the spinors U^{p,s) satisfy the Dirac equations — M±)U^{p) = 0, while the spinors U^{p,s) 
and U^^{p, s) satisfy the Rarita-Schwinger equations (fy— M±)Uj^{p) = 0 and (j/— M±)U^,,{p) = 0, 
and the relations l^"U^{p,s) = 0 , p^^U^{p,s) = 0 , l^U^Ap,s) = 0 , p^"U^^{p,s) = 0 , U^^{p,s) = 
U^^{p,s), respectively. On the other hand, the currents J^(0) and also couple potentially 

to the i , I"'" and ^ , I"*", | hidden-charm pentaquark states Pr, Pt and Pf, P^, Pr, 
respectively, 


(o|j^(o)|pr(p)) 

2 

(0|J^(0)|P^(p)) 

2 

(o|j;.fyo)|pr(p)) 

2 

(0|V(0)|P^(p)) 

2 

(o|v(o)|pr(p)) 

2 


fiPtiilbU {p,s), 

2 

>^t‘il5U^{p, s) , 

2 

9iP),Puil5U~{p,s) , 

2 

fp'l5[PtiU+{p,s)+p„U+{p,s)] , 

Xlij5U~^{p,s) , 


(19) 


( 20 ) 


the spinors (p, s) and P+ (p, s) s) and U^^,{p, s)) have analogous properties, and the pole 

residues Af ,5 ft ,3 7 ^ 0 and gt 7 ^ 0. 

2/2 2/2 2 

We insert a complete set of intermediate pentaquark states with the same quantum numbers as 
the current operators J^(x), i^ 5 J^{x), J^v{x) and i^^J^uix) into the correlation functions Ii^i,(p) 
and Xlfiuap (p) to obtain the hadronic representation (STJ [32] . After isolating the pole terms of the 
lowest states of the hidden-charm pentaquark states, we obtain the following results: 


n/xi/(p) 


Ilpi^a/S (p) 


A 


2 fl + M- 


2 Mf — 


—Pliy + 


+A 


_,_2 ^ — M+ 

I — p^ 


- 9 ^. 1 - 


l^Llu ‘^PfJ.Pv _ Pp,l,3 - Pul^ \ 

3 ) 

, lp.lv , 


pn 




Ml-p’^ 


PpPv + fl 


Ml - 


3p2 

■PpPv 


3v^ 


( 21 ) 


^ + 2 ^ + M+ 

I M2-p2 


7.7« + 


9pa9vP P 9pP9va 9 pv9olP f ( , IpPa laPp PpPa \ ~ 

- I ipia + - p= -::5- I 9v/B 


2 5 

IvPa - laPv PvPa 




9plS 


+ A 


_ 2 ^ — M- 

I Ml — p^ 


P 

9pot9vp P 9pP9voL 9 pv9olP 1 


1 / , IvPa - laPv PvPa 

--pr Ivla P 


P/. 

+/.+ 


10 

_ 2 P M_ 




5 

5/^/3 


10 


Ipla P 


7/.iPa - laPp PpPa 


\/p^ 


p^ 


9v/3 


Ml-p"^ 

2 ^ — M_|_ 
i Ml-p 2 

PM+ 
-2 



Ivlp 

‘ 2 pvP 0 

PvlP - pplv 

PpPa 1 gvP ~r 

3 

3p2 

3v^ 


Ivlfl 

, 2 p„p 0 

Pvlp - Pplv 

PfiPa 1 yiifi \ 

3 

H 

to 

3\/p2 


P • ■ 


P • 


, 2 |ipM_|_ _ 2 '^—M_ 

+5i Ml - p2 Pi^P''P<^PP + Ml - p2 Pt^P''P<-PP + • • • ’ 


( 22 ) 


where 5 ^ 1 / = the M± are the masses of the lowest pentaquark states with the parity P 

respectively, and the Af ,5, ft/s and gt are the corresponding pole residues. In calculations, we 

2/2 'll 1 2 
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have used the following summations [33) . 


s 


^ ^ , 7 m 7 ^ , Pi.lv-Pvl^ 

&( + Afa)(-!,,. + —+ - 




9iia9v0 H“ Pii/Spva. PfivPa^ 1 


7.7« + 


10 


71-7/3 


2 5 

IvPa IolPv PvPot I 

H p'" ) 

luPp - IpPu PuPp\ 


10 


7/.7q 


\/P 


5/^/3 I 7ai7/3 + 


5aio 


7/^P« - 7aPAt _ Pfi.Pa \ 

H P^ ) 

'JlJ.Pp - IpPfi PfiPp 


(23) 


Pup 


VI 


(24) 


and = M"^ on the mass-shell. 

We can rewrite the correlation functions II^j^ (p) and (p) into the following form according 

to Lorentz covariance, 


n^i/(p) = n|(p^) (-5^1,)-hn|(p2)7^7^-hn|(p^) (p^7,.-p^7^)-f Hi |(p2)p^pj,, (25) 

n^i-a /3 (p) = ns (p 2 ) gMagt -/3 + Pt^PPua + n| (p^) (p^^PaPp PaPPf,P,.) 

^ Z ^ 2 

+n| (p^) {PiialulP + PiiPlvla + Pvahilp + PuPl^ila) 

+n| (P^) [p,./3 (7MPa - laP^l) + Pua {liiPp - IpPfi) + Pm/ 3 (7i-Pa “ laPv) 

+Plia {ivPp - IPPv)] 

+n| J (p^) [p^iaPuPp + P,j,pPuPa + Pi^aPfiPp + PuPPiiPa) 

+n| J (p^) {np.laPvPp + ItilpPvPa + l^laPiiPp + IvlpP^iPa) 

+n| j(p^) [(7MPa - laPfi)PvPp + (7mP/3 “ lpPli)PyPa + (7i-Pa “ laPu) PiiPp 
+ {l>^PP -lPPv)PtiPa\ +Ilill{p‘^)Pfj.p„PaPp, (26) 

the subscripts i, | and | in the components 113 (p^), n|(p^), n|(p^), ni j(p^), n|(p^), n|(p^), 

n|(p 2 ), n|(p^), n|(p 2 ), n| 5 (p^), n| 5 (p^), n| ^(p^) and IIi 3 5 (p 2 ) denote the spins the pen- 
2 2 2 2’2 2’2 2’2 2 ’ 2'2 

taquark states, which means that the pentaquark states with J — and | have contributions. 

The components Hi 3 (p^), n| 5 (p^), n| 5 (p^), n| 5 (p^) and Hi 3 5 (p^) receive contributions 
2’2 2’2 2’2 2’2 2 ’ 2’2 

from more than one pentaquark state, so they can be neglected. We can rewrite 7 m 7 j- = P/iv — icTfiv^ 

then the components n|(p^), n|(p^), n|(p^) and n|(p^) are associated with tensor structures 
2 2 2 2 

which are antisymmetric in the Lorentz indexes p, a or fi. In calculations, we observe that 

such antisymmetric properties lead to smaller intervals of dimensions of the vacuum condensates, 

therefore worse QCD sum rules, so the components 11 3 (p^), n|(p^), n|(p^) and n|(p^) can also 

2 2 2 2 

be neglected. If we take the replacement Jfj,^(x) —>■ J^u{x) = J^u{x) — \Pi_iuJa‘^{x) to subtract 
the contributions of the J = \ pentaquark states, a lots of terms oc g^^, Pap disappear at the 
QCD side, and result in smaller intervals of dimensions of the vacuum condensates, so the com¬ 
ponents n 5 (p^) and n|(p^) are not the optimal choices to study the T = I pentaquark states. 

Now only the components n|(p^) and ns (p^) are left. The present conclusion is tentative, we can 

obtain definite conclusion by obtaining QCD sum rules based on the components n|(p^), n|(p^), 

2 2 

n 5 (p^), n|(p^), n|(p^) and n|(p^). in this article, we choose the tensor structures gau and 
2 2 2 2 

o± c ± 

PtiaPvp + PfipPva for analysis, thus separate the contributions of the | and | pentaquark states 
unambiguously, and tentatively assign the Pc(4380) and Pc(4450) to be the | and pentaquark 
states, respectively. 
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The current J^{x) has non-vanishing couplings with the scattering states pJ/ip^ A+Z)*°, pXci 
etc. In the following, we illustrate how to take into account the contributions of the intermediate 
baryon-meson loops to the correlation function n^i/(p), 




1 


(- M- - Sp(p) - (p) - Spxci (P) + ■ • ■ 

1 


\-2 


+*75 


T^-M+- S+j/^(p) - S|+ n.o(p) - ^txci ip) + 


-*75 -, (27) 


A+n*o' 


where the Af and M± are bare quantities to absorb the divergences in the self-energies Ep^^^(p), 
E^_,_^^g(p), Ep^__^(p), etc. The renormalized self-energies contribute a finite imaginary part to 
modify the dispersion relation, 




+ M_ ^_2 

M 2 +*v/^r_(p 2 ) t 


^-M+ _ 

M 2 +zy^r+(p 2 ) 



If we assign the Pc{4380) to be the </^ = | pentaquark state, the width r_(p2 = M^) = 
r p^( 438 o) = 205 ±18 + 86 MeV, which is much smaller than the width of the Zc(4200), r 2 ^( 42 oo) = 
370l7Q^™2 MeV. In Ref. [53], we observe that the finite width (even as large as 400 MeV) effect 
can be absorbed into the pole residue A 2 ^( 42 oo) safely, the intermediate meson-loops cannot affect 
the mass Mz^[i 200 ) significantly, so the zero width approximation in the hadronic spectral density 
works. The contributions of the intermediate baryon-meson loops to the correlation function 
II/i,ya/ 3 (p) can be studied analogously, furthermore, the width rp^( 445 o) is much smaller than the 
width rp^( 43 go). In this article, we take the zero width approximation, which will not impair the 
predictive ability significantly. 

Now we obtain the spectral densities at phenomenological side through the dispersion relation, 


Imlla (s) 
2 ^ ^ 

TT 


Imlls (s) 

2 ' ' 

TT 



m 2 ) 


M_A3^(5(s- 


m!) + M+At^^(s- 


Ml) - M+X+^S (s - 
Ml) - M_X-^^ 6 (s- 


Ml) , 

(29) 
Ml)] , 


= :^p|,p(s)+p|,p(s), 


(30) 


where the subscript H denotes the hadron side, then we introduce the weight function exp 
to obtain the QCD sum rules at the phenomenological side (or the hadron side). 


Uml 

/•so 

1 4m^ 

PSQ 

14m^ 

rSQ 

14m^ 


ds 

ds 

ds 

ds 


Vsplfjis) + p|^p(s) 
Vsp 3 p(s) - P 3 p(s) 
Vsplj^is) + p%uis) 
\fsp\ p(s) - p\ p(s) 


exp 

(-^) 

exp 

(-^) 

exp 

(-^) 

exp 

(-^) 


= 2 M-X 3 exp ( — 


Ml 

rp2 

, Ml 

= 2M+At exp ( -^ 

, Ml 

= 2M+At exp ( -^ 
, Ml 

= 2 M_A 5 exp ( 


{~w) 

(31) 

(32) 

(33) 

(34) 


where the sq are the continuum threshold parameters and the T'^ are the Borel parameters. We 
separate the contributions of the negative parity pentaquark states from that of the positive parity 
pentaquark states unambiguously. 
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In the following, we briefly outline the operator product expansion for the correlation functions 
(p) and Tlfivaisip) in perturbative QCD. We contract the u, d and c quark fields in the correlation 
functions and with Wick theorem, and obtain the results: 

{Tr [j 5 Dkk'{xh 5 CUf^,ix)C] Tr [ 7 ^C„„,(x) 7 .Ct/^™,(a:)C] CC^,,i-x)C 
-Tr [j5Dkk'ixh5CUl^,ix)Cj^Cnn'ixh.CUj'^,ix)C] CC^,,i-x)C} , (35) 


^fUfOC^ip) - 


^iia^ijk ^Imnl' a' j' k' m'n' 




{Tr [j 5 Dkk'{xh 5 CUjj,{x)C] Tr [y^Cnn'{xhaCU^^,{x)C] -f^CCj,^{-x)Cjp 
+Tr [y 5 Dkk'ix)-/ 5 CUj'j,{x)C] Tr ['y^Cnn'{xhaCU^^,{x)C] -t^CCl,^{-x)C'yp 
+Tr ["f^Dkk'{x)-t^CUjj,{x)C] Tr Yi^Cnn'{x)-^i3CU'^^,{x)C] j^CC^,^{-x)Cja 
+Tr [j5Dkk'{xh5CUjj,{x)C] Tr [y^Cnn'{x)'yi3CU^^,{x)C] 7^CC'J„(-a;)C'7„ 
-Tr [y5Dkk'ix)-f5CUZf{x)C-f^Cnn'{xhaCUj'^,{x)C] -i^CCl,^{-x)C-ip 


rnj 

tT 


jm' 

T 


-Tr [)zDkk'{x)-tzCUl,j,{x)C-tyCnn'{.x)-taCUj^,{x)C] Ji,CCj,^{-x)C-fp 
-Tr [j,Dkk'ixh5CUZ,>{x)Cj^C^n'{xhpCUf^,ix)C] ^.CCl,{-x)C^^ 

-Tr [iBDkk'{x)-i^CU'^j,{x)C-iyCnn'{x)-ii3CUj^,{x)C] 7 ^CC'^„(-a;)C' 7 „} , (36) 

where the Uij{x), Dij{x) and Cij(x) are the full it, d and c quark propagators respectively {Sij{x) = 


Sij{x) = 


Sij{qq) 5^jX^{qg^crCq) 


2tt'^x* 

-lidjc 


12 


192 


327r2a;2 


+ ■ 


(37) 


C^j{x) = 




(271)4 


d ke 


6ij PsG^pt'^j a°‘P{jt + rric) + (/^ + mc)(T°‘^ 


ji — rric 


{k? — mlY 




4 (A:^ — rrYY j 


(38) 


and t" = ^, the A" is the Gell-Mann matrix [35] , then compute the integrals both in the coordinate 
and momentum spaces to obtain the correlation functions 11^;^ (p) and Tlfj,^a^{p) therefore the QCD 
spectral densities plys QQjyis) and QQjjis) through the dispersion relation. In Eq.(37), we 

retain the term {gjUfiugi) comes from the Fierz re-arrangement of the (gidj) to absorb the gluons 
emitted from both the heavy quark lines and light quark lines to form {djgsG'^pt'^^a^uqi) so as to 
extract the mixed condensate {qggoGq). 

Once the analytical QCD spectral densities p|y 5 qcd^^) obtained, we 

can take the quark-hadron duality below the continuum thresholds sq and introduce the weight 










function exp to obtain the following QCD sum rules: 


2 M-X 3 ^ exp 
2 

2 M_|_At^ exp 
2 

2 M_|_At^ exp 
2 

2M_A7^ exp 



where 



'^p\,qcd(^) + P\,QCD 
'^P^,QCD^^) ~ P\,QCD 
'^P^,QCd(^) ~ P\,QCD 
'^P\,QCD^^) + P\,QCD 


exp 

:-^)' 

(39) 

exp 

:-^)- 

(40) 

exp 

:-^)' 

(41) 

exp 

:-^)' 

(42) 


P\,QCD^^) 

^|,QCn(®) 

P\,QGD^^) 

P\,QCDi^) 


Pqcd{s) ) 
2Pqcd('S) I 
'^cf^CDi.^) ! 
2™cPqcd('5) ) 


(43) 

(44) 


PQCd(s) = Po(s)+P3(s) + pl(s)+P5(s) + P6(s)+P8(s) + P9(s)+Pm(s): 

P°qcdW = pS(s)+^W + ^(s)+P^(s) + P^(s)+P°(s) + ^(s)+p°o(s), (45) 

the explicit expressions of the QCD spectral densities p\{s) and ^(s) with i = 0, 3, 4, 5, 6, 8, 9, 10 
are shown in the appendix. 

From Eqs.(39-44), we can see that if we set At = \/2At and Aj = v^A^, the four QCD sum 

2 2 2 2 

rules in Eqs.(39-42) are reduced to two QCD sum rules, the negative parity pentaquark states have 
degenerate masses, and the positive parity pentaquark states also have degenerate masses. The 
LHCb collaboration observe that the best fit leads to the spin-parity assignment (| , I"*") for the 
(Pc(4380), Pc(4450)), other assignments, such as (|^,f ) and (|''',| ), are also acceptable [15]. 
While Eqs.(39-44) indicate that the pentaquark states with the spin-parity (| , I”*") and (| , I"*") 
have degenerate masses, which contradicts with the assignments (I"*", | ) and (I"*", | ). 

In this article, we carry out the operator product expansion to the vacuum condensates up to 
dimension-10, and assume vacuum saturation for the higher dimension vacuum condensates, see 
Eqs.(35-38). We take the truncations n < 10 and fc < 1 in a consistent way, the operators of the 
orders 0(aj) with k> 1 are discarded. The condensates {g^GGG), { °'‘^^ ){sgsCrGs) have 

the dimensions 6, 8, 9 respectively, but they are the vacuum expectations of the operators of the 
order 0{as'^^), 0{ag), respectively. Furthermore, the numerical values of the condensates 

{qq){^GG) and (qq^i^GG) are very small, and accompanied by large denominators, and they 
are neglected safely. 

We differentiate Eqs.(39-42) with respect to w, then eliminate the pole residues Af.s. and 
obtain the QCD sum rules for the masses of the pentaquark states, 


Ml 


111,2 ds s [VspQcois) + Wc^cd('S)] exp (-^) 


111,2 ds yip^p^(s)-bme^p^(s) exp(-^) 


r^o 


4m2 


dss [\/spqcd(s) - exp (-^) 


IZ, ds 


VspqcDis)-mc'^gcois) exp(-^) 


(46) 

(47) 


9 


















where the M_ (M+) are the masses of the | (f*", f*") pentaquark states. Once the 

masses M± are obtained, we can take them as input parameters and obtain the pole residues from 

the QCD sum rules in Eqs.(39-42), the relations At = V^At and A 5 = \/2A7 hold. 

2 2 2 2 


3 Numerical results and discussions 


We take the vacuum condensates to be the standard values {qq) = —(0.24±0.01 GeV)^, {qggCrGq) = 
m1{qq), = (0.8±0.1) GeV^, = (0.33GeV)"^ at the energy scale g = I GeV [^[32]. The 

quark condensates and mixed quark condensates evolve with the renormalization group equation, 


as{Q) 

as(Ai) 


In the article, we take the 


( 99 )(a^) = ('?'?)(0) and {qgsaGq){g) = {qgsaGq){Q) 

MS mass mc{mc) = (1.275 ± 0.025) GeV from the Particle Data Group [34], and take into account 
the energy-scale dependence of the MS mass from the renormalization group equation, 


rricig) 

= mc{mc) 

’ a,ip) ' 
_as{mc). 

25 


as{g) 

1 

1 - 

61 log t bfilog^t- 

logt - 1 ) - 1 - 6062 ’ 

bot 

bl t 

btt^ \ 


(48) 


where t = log ^, 60 = 


S3 — 2nf 


h = 


153-19n/ 
2477 ^ 


bo = 


2857-^nf + ^7 


A = 213 MeV, 296 MeV 


12n ’ 247r2 > 128ir3 

and 339 MeV for the flavors n/ = 5, 4 and 3, respectively [34] . 

In Refs.[TH[lHI231l2SlllS], we study the acceptable energy scales of the QCD spectral densities 
for the hidden charm (bottom) tetraquark states and molecular (and molecule-like) states in the 


QCD sum rules in details for the first time, and suggest a formula g = ^^x/y/z ~ to 

determine the energy scales, where the V, V, Z denote the four-quark systems, and the Mg denotes 
the effective heavy quark masses. The effective mass Me = 1.8 GeV is the optimal value for the 
diquark-antidiquark type tetraquark states [HdilllS]. 

In this article, we use the diquark-diquark-antiquark model to construct the currents to inter¬ 
polate the hidden-charm pentaquark states, there also exists a cc quark pair. The hidden charm 
(or bottom) five-quark systems qqiq 2 QQ could be described by a double-well potential, just like the 
four-quark systems qq'QQ, see Eqs.(3-8) and related discussions in the introduction. The heavy 
five-quark states are also characterized by the effective heavy quark masses Mg and the virtuality 

V = Mp^ — (2Mg)2. The QCD sum rules have three typical energy scales /i^, T^, V^, we can 

also take the energy scale, = 0{T^) [T4l|26]. In this article, we can take the analogous 

formula. 


g = ^Ml^ - (2M,)2, (49) 

with the value Me = 1.8 GeV to determine the energy scales of the QGD spectral densities [I^[26] . 
and obtain the values g = 2.5 GeV and g = 2.6 GeV for the hidden charm pentaquark states 
Pc(4380) and Pc(4450), respectively. The energy scale formula can be rewritten as 

= (2Mc)^ -f ■ (50) 

In this article, we choose the Borel parameters and continuum threshold parameters Sq to 
satisfy the following criteria: 

1. Pole dominance at the phenomenological side; 

2 . Gonvergence of the operator product expansion; 

3. Appearance of the Borel platforms; 

4. Satisfying the energy scale formula. 
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Figure 1: The masses of the pentaquark states with variations of the threshold parameters Sq. 

In the QCD sum rules for the multiquark states, it is difficult to satisfy the criteria 1 and 
2. In previous work [Mj |25], we observed that the pole contributions can be taken as large 
as (50 — 70)% in the QCD sum rules for the diquark-antidiquark type tetraquark states qq'QQ 
{X,Y,Z), if the QCD spectral densities obey the energy scale formula /r = ^x/y/z ~ 

The operator product expansion converges more slowly in the QCD sum rules for the pentaquark 
states qqiq2QQ compared to that for the tetraquark states qq'QQ, so in this article, we choose 
smaller pole contributions, about (50 ± 10)%. For the tetraquark states qq'QQ [HI [25], the Borel 
platforms appear as the minimum values, and the platforms are very flat, but the Borel windows 
are small, = 0.4 GeV^, where the max and min denote the maximum and minimum 

values, respectively. For the three-quark baryons qq'Q, qQQ', QQ'Q" [30l|35], the Borel platforms 
do not appear as the minimum values, the predicted masses increase slowly with the increase of 
the Borel parameter, we determine the Borel windows by the criteria 1 and 2, the platforms are 
not very flat. In this article, we also choose small Borel windows = 0.4 GeV^, just 

like in the case of the tetraquark states, and obtain the platforms by requiring the uncertainties 
induced by the Borel parameters are about 1%. 

Now we search for the optimal Borel parameters and continuum threshold parameters sq 
according to the four criteria. The resulting Borel parameters, continuum threshold parameters, 
energy scales, pole contributions are shown explicitly in Table I. Furthermore, the contributions 
of the vacuum condensates of dimension 10 are less than 5%, the operator product expansion is 
convergent. So the four criteria of the QCD sum rules are satisfied, we expect to obtain reasonable 
predictions. From Table 1, we can see that the values ^so = Mp^^s- 5 +^ + (0.6 — 0.8) GeV ( or 

^ = (26 ± 1) GeV^, ^ = (27 ± 1) GeV^ ) can lead to satisfactory results. 

In Fig.l, we plot the predicted masses with variation of the threshold parameters sq, where 
we assign the Pc(4380) and Pc(4450) to be the | and I"*" pentaquark states, respectively. From 
the figure, we can see that the predicted masses increase slowly with (or are not sensitive to) the 
threshold parameters sq for central values of other parameters. 

In Refs. [301135], we study the = \ and | heavy, doubly-heavy and triply-heavy baryon 
states systematically with the QGD sum rules by subtracting the contributions from the corre¬ 
sponding and heavy, doubly-heavy and triply-heavy baryon states, the continuum 

threshold parameters yGo = Mgr-|-(0.6 —0.8) GeV work well, where subscript gr denotes the ground 
states. In the present case, the hidden charm pentaquark states carry a baryon number of one, i.e. 
they are doubly-heavy baryons. So the threshold parameters yGo = s+j + (0.6 — 0.8) GeV 


II 











r2(GeV^) 

^/i5(GeV) 

/i(GeV) 

pole 

Mp„(GeV) 

XpAGbY'^) 

m~) 

3.3- 3.7 

5.10 ±0.10 

2.5 

(40 - 61)% 

4.38 ±0.13 

(1.55 ±0.28) X 10-3 

^c(r) 

3.1-3.5 

5.15 ±0.10 

2.6 

(40 - 63)% 

4.44 ±0.14 

(0.84 ±0.17) X 10-3 


Table 1: The Bor el parameters, continuum threshold parameters, energy scales, pole contributions, 
masses and pole residues of the pentaquark states. 
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Figure 2: The masses of the pentaquark states with variations of the Borel parameters T^. 


make sense. One may worry that there exist some contaminations from the higher resonances, the 
upper bounds of the factors exp (—yS-) are about 0.0007 and 0.0004 in the QCD sum rules for 
the Pc (4380) and Pc (4450), respectively, if we take the largest values of the continuum threshold 
parameters, so the contaminations are greatly suppressed and can be neglected safely. 

We take into account all uncertainties of the input parameters, and obtain the values of the 
masses and pole residues of the | and | hidden-charm pentaquark states, which are shown 
in Figs.2-3 and Table 1. The QCD sum rules in Eqs.(39-42) and Eqs.(46-47) indicate that the 
pentaquark states with the spin-parity (| , |^) and (| , I"*") have degenerate masses, and At) = 

and A7 = \/2A7- Naively, we expect that additional one unit spin or P-wave can lead 
2 2 2 

to larger masses, so M 5 + > Ms-, while the relation M 3 + > Ms- needs detailed and refined 
2 2 2 2 

analysis to obtain the answer ”yes” or ”no”. It is sensible to assign the Pc(4380) and Pc(4450) 
to be the | and I"*" pentaquark states, respectively. However, the assignment (| , I"*") of the 
(Pc(4380), Pc(4450)) is not excluded. 

From Table 1, we can see that the present predictions Mp^( 438 o) = 4.38 ± 0.13 GeV and 
Mp,( 445 Q) = 4.44 ± 0.14 GeV are in good agreement with the experimental data of the LHCb 
collaboration, Mp^( 438 o) = 4380 ± 8 ± 29 MeV and Mp_,( 445 o) = 4449.8 ± 1.7 ± 2.5 MeV [15]. The 
present predictions support assigning the Pc(4380) and Pc(4450) to be the | and hidden charm 
pentaquark states, respectively, which are consistent with the assignments that the Pc(4380) and 
Pc(4450) are diquark-diquark-antiquark type pentaquark states [15] or the diquark-triquark type 
pentaquark states m- 

In this article, we take the energy scale formula /r = y^Mp^ — (2Mc)^ to determine the energy 
scales of the QCD spectral densities. The pole contributions are about (40 — 60)%, and the 
contributions of the vacuum condensates of dimension 10 are less than 5%, the two criteria (pole 
dominance at the phenomenological side and convergence of the operator product expansion) of 
the conventional QCD sum rules can be satisfied, so we expect to make reasonable predictions. In 
subsequent works, we extend the present work to study the ^ and | hidden-charm pentaquark 
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Figure 3: 


The pole residues of the pentaquark states with variations of the Borel parameters T^. 


states in a systematic way [36], where the energy scale formula /r = ~ (2Mc)^ serves as an 

additional constraint on the predicted masses. The typical energy scales, which characterize the 
five-quark systems qiq 2 q 3 cc and serve as the optimal energy scales of the QCD spectral densities, 
are not independent of the masses of the five-quark systems qiq 2 q?,cc. All the predictions can be 
confronted to the experimental data in the future. 

The diquark-diquark-antiquark type current with special quantum numbers couples potentially 
to special pentaquark states according to the tensor analysis in Eqs.(21-22) and Eqs.(25-26). The 
current can be re-arranged both in the color and Dirac-spinor spaces, and changed to a current 
as a special superposition of the color singlet baryon-meson type currents. The baryon-meson 
type currents couple potentially to the baryon-meson pairs. The diquark-diquark-antiquark type 
pentaquark state can be taken as a special superposition of a series of baryon-meson pairs, and 
embodies the net effects. The decays to its components (baryon-meson pairs) are Okubo-Zweig- 
lizuka super-allowed, but the re-arrangements in the color-space are non-trivial m- 

In the following, we perform Fierz re-arrangement to the currents and J^i, both in the color 
and Dirac-spinor spaces to obtain the results, 

= ^ Sccj ^ u + ^5757707^750- ^57^750077577- ^57^757707750 

-^57^0 077 - ^ 57^77 CO - ^5crAAiC 07^77 - ^Saxp,ucj^c+ ^Saxf,^5Ccj^^5U 

i 11 

+ -Saxfj.'f5UC'y^'f5C+ -Saxrjp.cca^'^u + -Saxr'fp.uca^'^c, (51) 

4 00 


Jfiiy — 


[guxiii + g^,xiu)cc-i^u -f + giixl^) u 07^0 


5 {guxl^l + gtixiv) 75007^7577 - -^5 {gyxl^^ + g^^xju) 75^07^750 




2x/2 


2x/2 


+ ^y^5 (7^crAr7i' + 7i^orAr77i) cccr'^'^77 -l- ^;^5 {■^fi.axTlu + 7i'0'at7ai) ^ ca^'^c , (52) 


where we take the replacement 


liu ^ 


J, 


/J.lv' 


J, 


1 


^ ^ilaAikAmn T ? 

= £ UjCj^dk 


1 ^ uC U^C"/^Cn"/^CCg "/ CCg 


(53) 
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to subtract the contribution of the spin-i pentaquark state, and use the notations 5rc = C^^djTck 

and iSFu = C^^djTuk for simplicity, here the F denotes the Dirac matrixes. 

The components 5(a;)Fc(a;)c(a;)F'u(a;) and 5(x)Fu(a;)c(x)F'c(a;) couple potentially to the baryon- 
meson pairs. The revelent thresholds are Mj/^p = 4.035 GeV, Mp^p = 3.922 GeV, M^^7v(i44o) = 

4.414 GeV, = 4.353 GeV, M^+qo = 4.151 GeV, M^+p^.o = 4.293 GeV, Mh,p = 4.463 GeV, 

^Xcip = 4.449GeV, and -Wa+( 2595 )do = 4.457GeV [33]. After taking into account the currents- 
hadrons duality, we obtain the Okubo-Zweig-lizuka super-allowed decays, 

P,(4380) ^ pJ/V-, A+D*° , P77, , A+.D° , px,o , (54) 

Pe(4450) ^ pJ/iA, A+D*° , P77, , A+D° , iV(1440)??e. (55) 

We can search for the Pc(4380) and Pc(4450) in the A+Zi*°, pr]c, A+Zi°, pxco, N{lM0)r]c mass 
distributions in the future, which may shed light on the nature of those pentaquark states. 


4 Conclusion 

In this article, we construct the diquark-diquark-antiquark type interpolating currents, and study 
the masses and pole residues of the | and | hidden-charm pentaquark states in details with the 
QGD sum rules by calculating the contributions of the vacuum condensates up to dimension-10 
in the operator product expansion. In calculations, we use the formula p = — (2Mc)^ to 

determine the energy scales of the QGD spectral densities. The present predictions favor assigning 
the Pc(4380) and Pc(4450) to be the | and I"*" pentaquark states, respectively. The pole residues 
can be taken as basic input parameters to study relevant processes of the pentaquark states with 
the three-point QGD sum rules. 
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Appendix 

The QGD spectral densities pl{s) and ^(s) with z = 0, 3, 4, 5, 6, 8, 9, 10 of the pentaquark states, 

^ 4915207r8 / - V - zf [s - mlf {7 s - 2ml) ^ 

Po(s) = J dydz{y + z){l-y-z)*{s-miy {6s-ml) , (56) 

P3(s) = -^^j^Jdydz{y + z){l-y-zf{s-mlf, 

^(s) = - J i^-y- (g -mlf , (57) 


14 







pUs) = 


= 


737287r6 
19 

'707788^ 

13 

'3932167r6 


J ^ (1 - y - {s - ml) ( 2 s - ml) 

^Os^) J (^y^^^l_y_ ^)3 _ —2^2 

^Os^G) J ^y^^y^(^i_y_ (s-m^)^ (5 s - 2 to^ ) , 


m; 


2949127r6 

1 

'2949127r6 

19 

' 11796487r6 
13 

' 7864327r6 


^ asGG ^ 


J (^ + ^ + ^ + ^) (1 - y - - ™c) (3s - ml) 

^Os^) J ^l_y_ 2)4 _ —2^2 _ _2^ 

^Os^) J {1-y- z)^ (s - mlf (2s - ml) 

J ^y^^ (y + -y- {s- ml)^ ( 4 s - ml), (58) 


pUs) = 


Mis) = 


mcjqgsCrGq) 
20487r6 
mcjqgsCrGq) 
655367r® 
mcjqgsCrGq) 
983047r6 ^ 

imcjqgsCrGq) 
327687r6 
mcjqgsfjGq) f 

10247r6 J 

mcjqgsCrGq) 
655367r6 J 
mcjqgsCrGq) 
983047r6 J 
imcjqgsCrGq) 
163847r6 


dydz jy + z)(l - y - z) {s - mlY 


dydz ( ~ + ~ ) (1 ~ y ~ (s — mlY 


z 

z y 

j dydz “ y “ (® “ ^c)" 

J dydz {y + z){l-y-z){s- ml)"^ , 
dydz jl — y — z) (^s — ml)"^ 

dydz {1-y- z){s- mlY , 


(59) 


Peis) = J dydz yz{l - y - z) [s - ml) ( 2 s - ml) 


jm)^ 


^(s) = J dydz jy + z){l - y - z) {s - ml) (3s - ml) 


(60) 


Psis) = - 


Mis) = - 


i^jqq)jqgscrGq) 

61447r4 

jqq)jqgscrGq) f 

122887r4 J 

s^jqq)jqgscrGq) 

122887r4 

jqq)jqgscrGq) f 
61447r4 / 


J dydz yz (3s — 2rril) 

dydz jy + z){l — y — z) (5s — 4:ml) , 

J dydz jy + z) ( 2 s — ml) 

dydz {1 — y — z) (4s — 3mg) , 


(61) 
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Vf 


pI{s) 


mc{qqf 

1447r2 



dy, 


mcjqq)^ 

727r2 



Plo('®) ~ 

19{qgsaGq)‘^ 

245767r4 J 

fVf 

1 dyy{l-y)[2 + mlS{s-ml)] 

Vi 


17 (qgsaGq)'^ 
4423687r4 

J dydz {y J- z) [4 J- J (s - ml)] 

= 

19{qgsaGqy 
491527r4 j 

fVf 

1 dy [l+mis {s- ml)] 

Vi 


17 (qgsaGq)'^ 
^ 2211847r4 

J dydz [3 + ml 6 [s — ml)] , 


(62) 


(63) 


where / dydz = dy ^ dz,yf 


= A’f:fdy^fodyJ. 


m" 


appear. 


i-y 


_ l + y^l-4mg/s _ l-v^l-4mg/g _ 2 _ fa+z)m) 

2 ’ W* 2 > ys-m'j ’ '"■£ yz ) 

di: — >■ ^ dz when the 6 functions <5 (s — and <5 (s — fh^) 
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